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ABSTRACT 


When  a  transversely  isotropic  elastic  body  that  contains  a  notch  or 
a  crack  is  under  an  axisynmetric  deformation,  it  is  shown  that  the 
eigenfunction  solution  near  the  singular  point  is  in  the  fora  of  a  power 


unction  s 

r*o  J -  *4-'-  <"• 

series  .  ,  (gpVf ,  (sp>r£$  ,  4^jf  *  Cvfr.i)  •  •  •  i 


in  which  {gif)  is  the 


fk»  vAo  .  bi¬ 

polar  coordinate  with  origin  at  the  singular  point  and  (&/  is  the  eigenva¬ 


lue,  or  the  order  of  singularity.  A  difficulty  arises  when 3  a*  well  as 


5+k  where  k  is  a  positive  integer  is  also  an  eigenvalue.  In  this  case 
the  higher  order  terns  of  the  series  solution  nay  not  exist.  A  nodified 
solution  is  required  and  is  presented  here.w^The  nodified  solution  has 
the  new  terns  p^+k  (lnp)  Ft  (^3)  ,  p^+k+1  (lnp)  F2  (^,  iT^Vv-^As  an  applica¬ 
tion,  we  consider  the  stresses  near  a  broken  fiber  in  a  composite  which 
is  under  an  axisynmetric  deformation.  The  interface  between  the  broken 
fiber  and  the  matrix  also  suffers  a  delanination.  This  creates  stress 
singularities  at  several  points  sone  of  which  require  the  nodified 
eigenfunctions  presented  here. 
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Chapter  I 
INTRODUCTION 

Even  though  light  weight,  high  atrength  compoaitea  have  been  widely 
uaed  in  the  induatriea,  rigoroua  enelyaia  of  atreaa  diatribution  in  e 
coapoaite  which  conteina  deleainetiona  end/or  broken  fibera  ere  atill 
lacking.  The  difficultiea  are  due  to  the  preaence  of  atreaa  aingulari- 
tiea  at  the  aingular  pointa  auch  aa  the  interface  crack  tip  and  the  edg- 
ea  of  the  broken  fibera.  Accurate  predictiona  of  the  atreaaea  near  theae 
aingular  pointa  are  iaiportant  not  only  for  atudiea  of  fracture  behavior 
of  materiala  but  alao  for  atudiea  of  general  atreaa  analyaia.  In  find¬ 
ing  atreaa  diatribution  in  an  entire  apeciaien  numerically  by  a  finite 
element  acheme,  one  may  uae  regular  finite  elementa  everywhere  except  at 
the  aingular  pointa.  At  theae  aingular  pointa,  apecial  elementa  are 
uaed  in  which  the  aingular  nature  of  the  atreaa  ia  given  by  an  analyt¬ 
ical  expreaaion. 

The  problem  of  finding  the  atreaa  a ingulari t iea  at  the  apex  of  an 
iaotropic  elaatic  wedge  or  notch  waa  firat  conaidered  by  Knein  (1926) 
and  Williama  (1952)  in  which  they  aaaume  that  the  atreaa  diatribution 
under  a  plane-atreaa  or  plane-atrain  deformation  can  be  expreaaed  in 
terma  of  a  aeriea  of  eigenfunct iona  of  the  form  p^f(^,5)  where  p  ia  the 
radial  diatance  from  the  apex  and  the  f  ia  a  function  of  the  polar  angle 
and  the  eigenvalue  5.  For  given  wedge  angle  and  homogeneoua  boundary 
conditiona  at  the  aidea  of  the  wedge,  there  are  in  general  infinitely 
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■any  eigenvalues  &  and  tha  associated  eigenfunctions  p^£(\p,5).  Particu¬ 
larly  iaportant  in  applications  is  when  one  or  aore  of  the  5's  is  neg¬ 
ative  and  the  stress  is  singular  at  the  apex.  It  was  shown  that  for  a 
two  diaensional  body  under  an  external  loading  the  negative  6  appears 
when  the  wedge  angle  is  larger  than  v  and  5  “  -1/2  is  a  double  root  when 
the  wedge  angle  is  2n  (i.e.  the  case  of  crack).  The  technique  is  applied 
to  a  crack  along  (Williaas  1959,  England  1965),  and  normal  (Zak  and  Wil- 
liaas  1963,  Cook  and  Erdogan  1972)  to  the  interface  and  to  other  geome¬ 
tries  of  isotropic  composites  (Bogy  1970,  1971,  Bogy  and  Wang  1971, 
Erdogan  and  Gupta  1972,  Delale  et  al  1984).  A  systematic  derivation  of 
the  equation  for  finding  the  singularity  3  was  given  by  Dempsey  and  Sin¬ 
clair  (1981). 

Investigation  of  associated  problems  for  anisotropic  materials  was 
started  by  Sih  et  al  (1965)  and  has  become  active  only  in  the  last  dec¬ 
ade  (for  example.  Bogy  1972,  Kuo  and  Bogy  1974,  Delale  and  Erdogan  1979, 
Sih  and  Chen  1981,  Hoenig  1982).  However,  these  studies  are  limited  to 
two-dimensional  singular  points.  There  are  singular  points  which  are 
three-dimensional.  Three-dimensional  singularity  analysis  for  isotropic 
materials  was  first  performed  analytically  by  Benthem  (1977,1980)  and 
Kawai  et  al  (1977)  and  numerically  by  Bazant  (1974)  and  Bazant  and 
Estenssoro  (1977,1980).  Extension  to  anisotropic  materials  and  compos¬ 
ites  was  considered  recently  by  Somaratna  and  Ting  (1986A.B).  In  Soma- 
ratna  and  Ting  (1986B)  finite  element  schemes  are  employed  to  determine 
the  order  of  singularity  at  a  three-dimensional  singular  point  of  any 
geometry.  In  the  other  paper  of  Somaratna  and  Ting  (1986A)  the  order  of 
singularity  is  determined  analytically  for  the  special  case  of  tran- 
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aversely  isotropic  materials  under  an  axisymmetric  deformation.  The  sin¬ 
gular  point  is  assumed  to  locate  on  the  axis  of  symmetry.  In  present 
study,  we  consider  the  ease  in  which  the  singular  point  is  not  on  the 
axis  of  symmetry. 

Fig. 1  shows  the  cross  section  of  an  axisymmetric  body  under  an  axi¬ 
symmetric  deformation.  The  material  is  assumed  to  be  transversely  iso¬ 
tropic  with  the  z-axis  being  the  axis  of  symmetry.  We  are  interested  in 
the  stresses  near  the  singular  point  R.  The  associated  problem  for  iso¬ 
tropic  materials  was  investigated  by  Delale  and  Erdogan  (1981)  and 
Delale  et  al  (1984) .  However,  their  objectives  are  different  from  ours 
and  hence  their  series  solution  is  different  from  the  one  presented 
here. 

A  series  solution  for  the  problem  is  developed  in  Chapter  II. 
After  presenting  the  basic  equations  for  transversely  isotropic  materi¬ 
als  under  an  axisymmetric  deformation  in  Section  2.1,  the  general  solu¬ 
tion  in  the  form  of  a  power  series  in  p  is  presented  in  Section  2.2. 
Application  of  the  stress-free  conditions  at  the  sides  of  the  wedge 
leads  to  equations  for  the  eigenvalue  6  and  the  coefficients  in  the  pow¬ 
er  series.  This  is  presented  in  Section  2.3.  It  is  seen  that  the  eigen¬ 
function  associated  with  an  eigenvalue  6  no  longer  contains  a  single 

« 5  e  /.I.  t  ,  .t.«  k..  »k.  ,  _ _  -5+ 1  i 


term  p° f(^,$).  It  also  has  the  terms  pd+1f r Ws 5) ,  p^+^ f2(^,5)  ... 
Therefore,  the  inclusion  of  the  second  and  higher-order  terms  in  the 
special  element  is  not  simply  the  inclusion  of  the  eigenfunctions  asso¬ 
ciated  with  the  subsequent  smallest  eigenvalue  <5.  A  similar  situation 
occurs  in  wedge  with  curved  sides  under  a  2-dimensional  deformation 
(Ting  1985).  The  derivations  presented  in  Section  2.2  and  2.3  are  for 


% 


A 


h 


I 


& 


,\ 


the  general  case  in  which  the  two  eigenvaluea  p,,  p2  of  the  elasticity 
constants  are  distinct.  The  degenerate  case  in  which  pt  “  p2  (of  which 
the  isotropic  material  is  a  special  case)  is  discussed  in  Section  2.4. 
For  2-dimensional  deformations,  the  displacement  of  the  singular  point  R 
(Fig.l)  can  be  ignored  for  the  singularity  analysis.  For  axisymmetric 
deformations,  one  cannot  ignore  the  displacement  of  the  singular  point  R 
in  the  r-direction.  A  particular  solution  associated  with  the  displace¬ 
ment  of  the  singular  point  is  presented  in  Section  2.5.  A  difficulty 
arises  when  8  as  well  as  5+k,  where  k  is  a  positive  integer,  is  an 
eigenvalue.  In  this  case  the  higher  order  terms  of  the  series  solution 
of  the  eigenfunction  cannot  always  be  determined.  A  modified  solution 
is  required  and  is  presented  in  Chapter  III.  We  can  see  in  Section  3.1 
and  Appendix  B  that  the  modified  eigenfunction  solution  has  the  new 
terms  (lnp)  Ft  (^,  5)  ,  p<5+k+l  (inp)F?  (^,  5)  ...  Application  of  the 
stress-free  boundary  condition  is  presented  in  Section  3.2.  The  solu¬ 
tions  of  Chapter  II  are  then  applied  to  composite  materials  in  Chapter 
IV.  The  equations  for  general  transversely  isotropic  materials  and 
degenerate  materials  are  presented  in  Section  4.2  and  4.3,  respectively, 
and  numerical  examples  are  given  in  Section  4.4.  For  the  singular  point 
which  is  the  tip  of  an  interface  crack  the  displacement  was  found  to  be 
oscillatory.  This  implies  that  the  two  crack  surfaces  inter-penetrate 
each  other.  To  avoid  the  unrealistic  phenomenon  a  contact  zone  near  the 
crack  tip  is  introduced  in  Chapter  V.  In  Chapter  VI  the  formulas  for 
singularities  at  an  interface  crack  with  contact  zone  are  derived  by 
using  the  Stroh  formalism.  This  alternate  approach  offers  an  analytical 
solution  for  the  singularity  8  which  agrees  with  the  numerical  results 
obtained  in  Chapter  V. 


EIGENFUNCTION  FOR  AXI-SYMMETRIC  DEFORMATIONS 


2.1  MATHEMATICAL  FORHJLATION 

Let  (r,0,z)  be  a  cylindrical  coordinate  system  with  the  z-axis  as 
the  axis  of  material  symmetry  and  let  (ur,u£,uz)  be  the  corresponding 
displacement  components.  He  assume  that  the  deformation  is  axisymmetric 
and  UQm  0  so  that  ur  and  uz  are  functions  of  r  and  z  only.  Introducing 
the  displacement  potential  4>( r,z)  which  gives  ur  and  uz  by  (Elliott 
1948,  Green  and  Zerna  1954  and  Kassir  and  Sih  1975) 


u 


r 


<>r 


» 


(2.1) 


where  m  is  a  constant  to  be  determined,  the  stresses  are  obtained  as 


b*<t>  dd>  d’4 

a_  ■  c  -  +  c  -  +  c  m  - 

r  )ri  r»r  13  bz* 


*<t>  a*d» 

O  0  “  C .  _  T—  +  C . .  —  +  c .  _  m 


(2.2) 


i  j  &ra  i 1 


rdr 


1  3 


*<t> 

O  m  C  —  +  c  —  +  c 
Z  13  ara  13  rfcr  33 


d3d> 

arz“  c«/1+m^  nr  » 

**  44  dr®z 


8z3 


in  which  c^j  are  the  elasticity  constants  for  the  transversely  isotropic 
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11*  ♦  _2*  -  i_  11*  0 

3r*  r^r  p*  3z* 


(2.3) 


where 


e,,  ♦  (l+m)i 


»c,,  +  <1+«)c44 


(2.4a) 


or,  equivalently. 


C1»  1  £l 

<C1»  +  C44>  p* 


C..  +  c. 


C  4  4  +  C,,  P* 


(2.4b) 


The  second  equality  of  (2.4a)  and  (2.4b)  respectively  yield 


-  2  [* 


c.  .C..  -  c 


— LA—  -1  1  a  ♦  1  -  0  , 

+  cA.)  J 


2c  (c  +  c 

44^1  1  4  4 


4  +  2  f  c,.L~n  2ci.tli.i . ]  D*  +  Lli 

L  2c.  e  J  P  e__ 


(2.5a) 


(2.5b) 


It  can  be  shorn  (Eshelby  et  al  1953)  that  p  cannot  be  real  if  the  strain 
energy  is  positive  definite.  Therefore  we  have  two  pairs  of  complex 
conjugates  for  p  which  will  be  denoted  by  p,,  p,,  p2  and  p2  where  an 
overbar  indicates  the  complex  conjugate.  The  associated  values  of  a  are 
denoted  by  m,t  a,,  m2,  m?  respectively.  From  (2.5a)  we  note  that 


-  1  . 


(2.5c) 


Since  (2.5b)  is  a  quadratic  equation  in  p *  with  real  coef f icient* ,  if  p, 
is  purely  imaginary  so  is  p?.  Then  m,,  u}  are  real  and  m,  -  m,,  m?  - 
If  p,  and  p2  are  not  purely  imaginary  we  may  choose 


p  ■  u  +  iv  ■  -o.  . 


o.  "  -u  +  iv  ■  -o 


where  u,  v  are  real.  In  this  case  a,  and  a2  are  complex  and  m,  “  m?.  In 
view  of  the  fact  that  the  equations  are  linear,  the  general  solution  for 
4  is  obtained  by  superposing  4>'s  associated  with  pt,  p , ,  pJt  p;.  We  will 
assuae  that  p,  *  p?.  The  degenerate  case  in  which  p,  “  p2  will  be  dis¬ 
cussed  in  Section  2.4  . 


2.2  EIGENFUNCTIONS  FOR  SMALL  ft 

Let  (r,z)*(e,0)  be  a  singular  point  which  may  be  the  apex  of  a 
wedge,  notch,  crack, or  the  tip  of  an  interface  crack.  We  now  consider 
the  case  in  which  a  **  0.  The  case  in  which  a  ■  0  has  been  studied  by 
Soaaratna  and  Ting  (1986a).  Using  the  singular  point  as  the  origin,  we 
def ine 


x  -  r-a  *  p  cout  ,  z  •  p  sin^  . 


(2.7) 


To  find  the  eigenfunction  for  4  that  is  valid  for  small  p,  we  rewrite 


equation  (2.3)  as 

a  l  a  *<b 

'  p7  Sz7  ' 


l  a 4> 
a+x  ax 


-  Y  2  (  -  )\ 
a  dx  s*0  a 


(2.8) 


.  *«»_  i  *<1).  L  *(«. 


<t>  m  <t> 


■  -io  (4)Vk)' 


*<k)-  , 


Z  -  x  ♦  pz  , 


(2.9a) 


(2.9b) 


(2.10) 


i  ^  ^  * 


i  »  «  «  t  4  »  I.*  e.» 
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where  2  is  the  eigenvalue  end  Ap ^  are  constants  to  be  determined.  Using 
equation  (2.7)  we  have 


xtz«+k-t+2.  pd+k+2 (cog^) t  j-S+k-t+2  f 


(2.11) 


f  “  cos^  +  p  sin^  . 


(2.12) 


Therefore  is  of  order  p^+*+Z.  By  substituting  equations  (2.9)  into 
(2.8)  and  equating  the  coefficients  of  xtz^+k_t+^,  £t  can  be  shown  that 
(see  Appendix  A) 


s(k)-  ^  ■  o  >  o)  • 


(2.13a) 


k(k). 


2t 


i  K(k'1)-  (tM)Aftu  . 


(t  -  k-1 . 1) 

Hence  the  only  unknowns  are  A^  (k*0,l,2,  ...)  and  5  which  will  be 

determined  from  the  boundary  conditions. 

Me  will  let  the  solution  given  by  equations  (2. 9)- (2. 13)  apply  to  p 
“  p,.  For  p  ■  p2,  p ,  and  p2  we  will  use  the  same  expressions  except  that 
A^  is  replaced  by  ,  C ^  and  ,  respectively.  Thus  the  general 
solution  for  i* 


*(k)«  S  {A^x^f’* 


5+k-t+2+B  (k)  vt75+k-t+24.r  (k)  „t;5+k-t+2 

+Dt(k)xtzf+k_t+2}  . 


(2. 14a) 


Z,  •  x  +  p,2  (*”1.2) 


(2.14b) 


Subtituting  (2.14)  into  (2.1),  we  obtain 


V*.  13) 


u*m  lt-oA* '  (i+k't+2)  ♦tJ0At+i(t+1)i  *tz}+k_t+1 

+  {  2  Bt<k>  (i+k-t+2)  ♦Vb^U+D)  .tjfi+k-t+i 
t-0  1  t-0  1  * 

+  {  2  C<k> (i+k-t+2)  ♦Vcftht+l))  xl2f+k-t+1 

t-0  *  t“0  1 

+  {J0Dt(k)  (a+k-t+2)  +k^Dt(k}(t+i)}  xt2|+k"t+1 


To  siaplify  th«  expression,  wo  will  lot 
A^O-  B^k)-  C^k)«  D^k>«  0,  if  t>k  , 
and  rowrito  equation  (2.15)  as 


u<k)-  t20  [At(k)  (a+k-t-2)  ♦  A^U+l)] 


ra+k-t+i 


(2.16) 


(2.17a) 


whoro  the  dots  stand  for  the  siailar  expressions  associated  with  p  -  pJt 
p,  and  p2.  Likewise,  we  have 


t20At(k)B,p,(a+k-t+2)  xtz}+k“t+1  ♦ 


(2.17b) 


In  substituting  equation  (2.14)  into  (2.2)  for  the  stresses,  we  first 
replace  the  terns  d4/(rdr)  in  (2.2)  by  da4/dx*  and  da4/&z*  using  equa¬ 
tion  (2.3).  He  then  have 


C^k)-  troI-At(k)c44(l+B,)pa(a+k-t  +  2)  (6+k-t  +  l) 
+At(^2(c,  ,-c,j)  (t+1)  (i+k-t+l) 

+At(ki(c,  ,-c,,)  (t+2)  (t^DJxtzf^-t 


(2.18a) 


♦  . . . 


+  ••• 


(2.18c) 


(c ,  ,-c , ,)  (t+1)  (5+k-t+l) 
-A$$(t+2)  (t+Dlxtzf'*"1  +  ... 


t,5+k-t 


a*k)“  t-o  *Atk)c44(1+*i>(a+k-t+2>(a+k't+1)xtZl 
oQ -tXQ  [A^kb+k-t+2)  +a£J$  (t+1)  3  c4 4  (1+m, )  p ,  (5+k-t+l)  xtzf+k't+ . .  P ' 18d) 


In  (2.18a)  and  (2.18c),  the  following  identities  which  are  obtained  from 
equation  (2.4a)  have  been  used: 


cii  +  -c44(i+«V  . 

cn  +  ca»*P*“  -c44(l+«)  • 


(2.19a) 

(2.19b) 


2.3  DETERMINATION  OF  6  AND  ... 

The  problem  reduces  to  the  determination  of  5  and  A^k)  ...  Dp^ . 
The  stress-free  boundary  conditions  it  ^  #  and  a'  are 


ffrsin^  -  <7rzcos^  -  0  ,  (2.20a) 

arz sin^  -  <rzcos^  •  0  . 


Written  in  matrix  notation,  we  have 

i 

JW  2  -  o  ,  (*  -  a,  or’)  ,  (2.20b) 

i 

B  where 


From  equations  (2.18)  we  may  write  a 


as,  using  equation  (2.11), 


(k) 


.Ck).  /*tt20(StW,itk)at(k)*TtW,J.k)3t(k)1»ut«.«.k)St(^) , 


(2.23) 


in  which 


St 


(k)  . 


A<k> 

C^k) 

D<k) 


(2.24) 


St(^,5)  -  c44 (5-t+2) (S-t+1) 


(l+m,)pj  * 

(l+m,)p,  * 

(1+m,)  * 


*  * 
*  * 
*  * 


QtGM) 


(2.25a) 


TtGM)  -  (t+D(«-t+l) 


2(c,rc12)  *  * 

c«4<1+»i>Pi  *  * 

0  .00 


Ut0M)  -  (t+2)(t+l)(c11-c1J) 


1  1 
0  0 
0  0 


* 

* 

0 


Qt  W.5) . 


1  1 
0  0 
0  0 


£t  W.S) 


(2.25b) 


(2.25c) 


In  equations  (2.25),  the  second  column  of  the  matrix  is  obtained  from 
the  first  column  by  replacing  p,  by  p2  (and  of  course  m,  by  m;)  .  The 
third  and  fourth  columns  are,  respectively,  the  complex  conjugate  of  the 
first  and  second  columns.  Qt(^, 6)  is  a  diagonal  matrix  given  by 
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Qt(^M)  -  (cos^diag  [  f*-t,  fj'1,  ( f 5 - c ,  (fj)0-t  ]  .  (2.26) 


<1^  of  equation  (2.24)  ia  related  to  g£k-^  by  using  equatons  (2.13). 
Notice  that  equation  (2.13a)  can  be  regarded  as  a  special  case  of  equa¬ 
tion  (2.13b)  if  we  use  equation  (2.16)  and  allow  t  ■  k  in  (2.13b).  Thus 


we  have,  for  kkl. 


3tk>"  s&‘)+  Kk-»-  <tn,,c;i], 


(2.27a) 


(t  -  k,k-l . 1) 


q** ' *  0  ,  i 


if  t>k  . 


(2.27b) 


As  in  (2.13)  the  only  unknowns  are  6  and  (  k  -  0,  1,  2  ...  ). 


Before  we  substitute  (2.21)-(2.23)  into  (2.20b),  we  rewrite  (2.23) 


as.  Baking  use  of  (2.27b) 


»(k)-  -**k  fs.W.WOsJ 


(2.28a) 


+t-l  ^5t  5+k)  +It- 1  a+k>  +Ht-2  6+k)  3  St  * 


where  we  have  defined 


Uj.(^,5+k)  ■  0  ,  if  t<0 


(2.28b) 


Now  substitution  of  (2.21),  (2.22),  and  (2.28a)  into  (2.20b)  for  ^  -  a 


and  a'  yields  the  following  equations  for  q^k)  : 


K(5)  qo(0)-  0  , 


(2.29a) 


K(a+k)S^k)-  -t21Wt(5+k)gt(k)  (k>l)  , 


(2.29b) 


For  •  nontrivial  aolution  of  g£k^ ,  wa  tat  fro*  (2.29a)  that 

|  K(«)  |  -  0  .  (2.31) 

Thua  5  ia  tha  aiganvalua  of  tha  matrix  K  and  is  tha  associatad 

aiganvactor.  With  &  and  g^  obtainad  from  (2.31)  and  (2.29a),  (2.29b) 
providas  g(k)  for  k  i  1  and  (2.27a)  givas  g^  for  1  5  t  5  k. 

Whan  i  ia  a  siapla  root  of  (2.31),  g^^  obtainad  from  (2.29)  is 
uniqua  up  to  an  arbitrary  aultiplicativa  constant,  g^  for  k  >  1 

obtainad  froa  (2.27a)  is  uniqua  in  taras  of  g^  providad  £+k  ia  not  a 
root  of  (2.31).  Tharafora,  whan  5  is  a  siapla  root  and  &+k  is  not  a  root 
of  (2.31),  tha  aiganfunction  4>  associatad  with  5  is  uniqua  up  to  an 
arbitrary  auliplicativa  constant.  If  S  is  a  aultipla  root,  say  a  doubla 
root  of  (2.31),  and  (2.29a)  providas  two  indapandant  g^^ »  wa  would  hava 
two  indapandant  aiganfunct ions  aach  of  which  is  uniqua  up  to  a  aultipli¬ 
cativa  constant  providad  d+k  is  not  a  root  of  (2.31).  Whan  5+k  is  also 
an  aiganvalua  of  K,  wa  saa  froa  (2.29b)  that  a  solution  for  q^k)  axists 
if  and  only  if 

iTJi  St«*»st<k)  ■  o 


(2.32a) 


(Hilderbrand  1954)  where  the  superscript  T  denotes  the  transpose  end 
is  the  left  eigenvector  of  K(3+k) 


LT  K(3+k)  -  0  . 


(2.32b) 


If  (2.32a)  holds,  g^k^  exists  but  is  not  unique.  However,  the  nonunique 
portion  of  g^  c*n  ilnor*d  because  that  portion  is  represented  by 
the  eigenfunction  associated  with  the  eigenvalue  3+ k.  An  exaaple  of  this 
case  in  a  related  problea  can  be  found  in  Deapsey  (1981) ,  Zwiers  et  al 
(1982),  and  Ting  and  Chou  (1985). 

If  3+k  is  an  eigenvalue  of  K  and  (2.32a)  does  not  hold,  a  solution 
for  does  not  exist.  In  this  case,  the  solution  for  4^)  cannot  be 
given  by  equation  (2.9b).  Instead,  we  use  the  folowing  aodified  solu¬ 
tion: 


*oo.  L  i  Aik)xtza+k't'‘'2  (2,33) 

9  *3  t-o  1 

in  which  is  now  assuaed  to  depend  on  3.  This  case  will  be  discussed 
in  Chapter  III.  Equation  (2.33)  can  also  be  used  for  second  independent 
solution  when  3  is  a  double  root  of  (2.31)  but  (2,29a)  provides  only  one 
independent  g^  , 

We  see  froa  equations  (2.22)  and  (2.28)  that  for  each  eigenvalue  3, 
the  stress  has  the  teras  pi+1ft(^,3),  pi+2f?(^,3),  ...  Thus 
the  eigenfunction  associated  with  an  eigenvalue  has  infinite  teras  for 
axisyaaetric  deformations.  We  also  see  thet  if  Re(3)  <  0,  the  stress  is 
singular  at  p  ■  0.  Thus  Re(3)  provides  the  order  of  singularity. 


ft- 


2.4  DEGENERATE  CASE  m,  -  id,  -  1 

When  p,  "  p ? ,  p  aust  be  purely  iaaginary.  .  This  follows  froa  equa¬ 
tion  (2.5b)  and  the  fact  that  p  cannot  be  real.  By  equation  (2.4b)  and 
(2.5c),  we  have  a,  *  a2  ■  1.  We  cannot  have  -1  because  this 
would  sake  p  real.  By  setting  a  •  1,  the  second  equality  of  equation 
(2.4a)  yields 


<«..♦  *c44)’  - 


eiicss 


(2.34) 


Hence  when  equation  (2.34)  is  satisfied,  p,  ”  p,  and  we  have  a  degener¬ 
ate  case.  The  five  aaterial  constants  are  now  reduced  to  four  by  equa¬ 
tion  (2.34).  Introducing  the  new  aaterial  constants  v,  m,  y,  and  3,  we 
let 


c„-  (  X  ♦  2m  )  3'. 
c jj*  (  X  ♦  2m  )/3\ 


c 


44 


M  , 


C 


1  » 


X  , 


c,,*  C,,-  2-VM  , 


in  which 


X  -  2m*'/ (  1  -  2v  )  . 


(2.35a) 

(2.35b) 

(2.35c) 

(2.35d) 

(2.35e) 

(2.35f) 


Equations  (2.35)  satisfy  (2.34).  With  (2.35),  equations  (2.5)  give 
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For  isotropic  materials  w«  have  6*7*1  and  v  and  m  are  the  Poisson's 
ratio  and  shear  aodulus,  respectively. 

In  a  degenerate  case  p,  *  pJt  the  terms  associated  with  B^  and 


are  identical,  respectively,  to  the  terms  associated  with  A ^  and 


Cfk>.  He  therefore  need  a  new  solution  for  B ^  and  D^K  This  can  be 
accomplished  by  replacing  the  coefficients  of  B^k^  and  D^k^  by  their 
derivatives  with  respect  to  p,  and  pa  (Ting  and  Chou  1981B,  Ting  1982). 
Thus,  for  instance,  equation  (2.14a)  becomes 


*<k>-  x  t(k)*t£5+k-t+2) 

t-0  1  c 

+  X  {Bik^zxtZi+k't+1+Dp^zxtZ5+k'"t+1}  (6+k-t+2) 
t-0  1  1 


(2.37) 


where,  since  pt  -  p},  we  have  omitted  the  subscripts  1  and  2  for  Z  and 
I.  Similarly,  equations  (2.17)  and  (2.18)  are  replaced  by  (noting  that 
m  in  the  B £k^  and  D^k)  terms  must  also  be  differentiated  with  respect  to 
p  by  using  equation  (2.4b)), 


urk^ “  t£0 [Atk\$+k-t+2) +A^k{ (t+1) ] xtZ5+k_t+1 

+  X  [B»k^4+k-t+2)+Bfk{ (t+1)] (i+k-t+l)zxtZ5+k-t 
t-0  1  1 1 


(2.38a) 


uzk>"  t5olAtk)0i*Btk) (3-4v)] (5+k-t+2)xtZ 

k 


a+k-t+i 


(2.38b) 


♦  X  B*(k)0i(a+k-t+2)  (a+k-t+Dzx^5*11-6  + 
t-0  1 


i 


3 


Hi 


A  V- 


-iff <k)-  2  f(A^k)e*-B^k)2v3i) (6+k-t+2) (6+k-t+l)  (2.39a) 

2m  t-0  1  l 

+A^2T(t+l)  (6+k-t+l) +A^y(t+2)  (t+l)]xtZi+k_t 

)t 

+  t20[Bt(k)e*  (6+k-t+2)  (5+k-t+l)+B£k^2y(t+l)  (6+k-t+l) 

+B^ir(t+2)  (t+i)]  (a+k-Ozxtz^k-1-1*  . . . 

-i  <rjk)-  2  {[At(k)  (e*-Y)-Bt(k)2v3i]  (5+k-t+2)  (6+k-t+l)  (2.39b) 

2m  t-0  u 

-A^2Y(t+l)  (6+k-t+l) -A^^(t+2)  (t+l)}«t2,+k”t 

♦  Z  [B^kte*-y)  (6+k-t+2)  (6+k-t+l)-B^k^27(t+l)  (6+k-t+l) 
-Bt(£jy(t+2)  (t+1)]  (6+k-t)zxtZ5+k‘t_1+  . . . 

_i  <r<k)-  I  - [A^k)  +B5k) 2 (1-v) i3-1]  (6+k-t+2)  (6+k-t+l) xtza+k"t  (2*39c) 
2m  z  t-0  c  1 

■t50Btk^6+k-t+2) (6+k-t+l) (6+k-t)zxtZ5+k_t“1+  . . . 

-1  a^k)-  2  {[A^k)0i-B^k) (1-2*)] (6+k-t+2)  (2.39d) 

2m  z  t-0  1  t 

+  [Af5)3i-Bt^(l-2v)]  (t+1)}  (6+k-t+l) xtZ*+k"t 
+^2^ [B^k^6+k-t+2)+B^kj (t+1)] (6+k-t+l) (6+k-t)3izxtZ5+k-t_1+  .. 

In  equations  (2.38)  and  (2.39),  the  dots  stand  for  the  C^k)  and  D^k) 
terms  that  are  obtained  from  the  A^k^  and  B^k^  terms  by  replacing  3i  and 
Z  by  ~3i  and  Z,  respectively.  Equations  (2.20)~(2.32)  remain  valid 


except  (2.25),  which  are  replaced  by 


stGM) 


-  2/u(6-t+2)  (5-t  +  l)  3i 


(2.40a) 


-2«'3i+3*  (6-t)f~^sin^  *  * 

-(l-2v)+(6-t)3if~^sin^  *  *  Qt(^,6), 

-2(l-v)i3_1-(6-t)r1fin^  *  * 
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Tt(^,5) 


-  2M(t+l) (5-t+l)  0i  - 


2y(8-t) f-*sin^  * 

(l-2v)+ (6-t)0if~^«in^  * 


(2.40b) 


ut(iM) 


|"l  (8-t) f'^sin^ 


(2.40c) 


-  2MY(t+2) (t+1)  0 


0  0  Qt  (*,  8) 


0  0 


Qt(^,8)  is  obtained  from  (2.26)  with  ft  ■  {2  m  f.  The  third  and  fourth 


coluans  of  the  aatrices  in  (2.40)  are  obtained  froa  the  first  and  second 


coluans  by  replacing  0i  and  f  by  -0i  and  ?,  respectively.  Equatii 


(2.27)  reaains  valid  because  the  order  of  the  differentiation  with 


respect  to  p  and  x  or  z  can  be  interchanged. 


2.5  PARTICULAR  SOLUTION  FOR  THE  DISPLACEMENT  OF  THE  S INSULAR  POINT 


For  a  two-diaensional  problea,  the  displaceaent  of  a  singular  point 


can  be  ignored  for  the  singularity  analysis.  For  an  axisyaaetric  defor- 


aation,  one  cannot  ignore  the  displaceaent  ur  of  a  singular  point.  We 


therefore  consider  (ur,uz^  ”  (u0,0)  at  the  singular  point  (r,z)  “  (a,0) 


where  u.  is  a  constant.  A  particular  solution  that  yields  this  displace¬ 


aent  is 


ur-  u0r/a, 


uz"  0, 


(2.40a) 


rr-  eg-  u0(c, ,+  c, ,) /a, 


azm  2U0C13/#* 


atzm  °- 


(2.40b) 


y 


To  satisfy  the  stress-free  boundary  conditions,  equation  (2.20b),  we 
superinpose  equation  (2.42)  to  (2.23)  with  6-0  and  write  the  stress  as 


z  “  Zo *  J0  px  Jo  {5t^.k)atk)+lt^*k)at+l+Ht^»k)at+i  }* 


Equation  (2.20b)  now  provides  the  following  equations  for  q 


(k) 


K(k)  q0(k)-  b(k>,  (k  -  0,  1,  2,  ...  ), 


in  which 


N(a) 

N(a’) 

•>* 


Zo 


t 


H(k)- -,5i  “t<k>stk).  <ki». 

and  (k)  is  defined  in  equation  (2.30b).  Equation  (2.44)  has 
solution  for  g^k^  if  k  is  not  an  eigenvalue  of  K  .  If  k  is, 
discussion  presented  near  the  end  of  Section  2.3  applies  here. 


(2.43) 


(2.44) 


(2.45a) 

(2.45b) 

a  unique 
then  the 
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Chapter  III 

MODIFIED  EIGENFUNCTIONS 


3.1  MODIFIED  SOWTIOM 

Whan  3  is  a  root  of  (2.29a),  lat  v  be  tha  saallast  positive  integer 
for  which  3+ij  ia  also  a  root  of  (2.29a).  Equation  (2.29b)  for  k  ■  ij  is 


-  -Jj  Jt  («*■>) s,(,>- 


(3.1) 


This  has  a  solution  for  unless 


iT,ii  '  »  . 


(3.2) 


where  L  is  the  left  eigenvector 


LT  K (3+n)  -  0  . 


(3.3) 


If  (3.2)  holds,  a  solution  for  go’’)  does  not  exist  and  the  expansion  for 
4  given  by  (2.9)  is  not  valid. 

To  obtain  a  valid  expansion  when  (3.2)  holds,  we  notice  that  if 
given  by  (2.9)  satisfies  (2.8)  so  does  &4/d£.  Therefore,  in  place  of 
(2.14a)  we  use 


*Ck)-  III  Z  A,00.'*4*1''*2*  ...  1 

00  t“  o  *■ 


(3.4) 


S 


-1 


.1 


. . .  .  .  .  . . j 


21 


ti  •*(  SI  SI  St 


in  which  A^)  ...  D^)  are  now  regarded  as  functions  of  b  (Zwiers  et  al 
1982).  If  we  carry  out  the  differentiation  in  (3. A),  we  will  have  terns 
of  the  fora  xaZ^  as  well  as  xaZ^(lnZ).  By  substituting  (3. A)  into  (2.9a) 
and  then  into  (2.8),  the  coefficients  of  xaZ^  and  xaZ^(lnZ)  oust  vanish. 
The  latter  leads  to  (2.27).  The  foraer  leads  to  the  following  equations 
which  can  also  be  obtained  by  differentiating  (2.27)  with  respect  to  6: 


St <k)"  Hr  st'-i'0*  2(i,,_t.2)  [  tst'(k"1)-<t*»s(ii)  1 
■  irrrnr  [  t3t<k"1)"  (t+,)s&  1  • 


(3.5a) 


for  t  “  k  ,k~l ,  ...  ,1  , 


(k)  . 


0  ,  for  t  >  k  , 


(3.5b) 


where  the  priae  denotes  differentiation  with  respect  to  6.  Hence  the 


ily  unknowns  are  and  q"  ,  (k  *»  0,  1,  2  ...) 


When  (3.2)  holds  and  v  is  the  smallest  positive  integer  for  which 


b+v  is  also  a  root  of  (2.29a),  we  may  choose 


*  0  ,  for  k  <  v  . 


(3.6) 


Substituting  (3. A)  into  (2.1)  and  (2.2)  and  carrying  out  the  differenti¬ 
ation,  we  obtain  new  expressions  for  the  displacements  and  stresses. 
This  is  presented  in  Appendix  B. 

for  the  degenerate  case  in  which  pf  “  p;,  we  use  instead  of  (3. A) 


*<k)-  1  (tJoUt<k)»t2i*k"t*2»Ct(k)KtZi*k't*2) 

.  J0[Bt<k)z»tZ{*k't*1.Dt(k)z.tii*k't‘I]  (5+k- 1*2) i 


(3.7) 


.cot 


k*  -“'i  w'e to'*  .Ti  / «  'w  *.  .  k  .  -  .*•  .  «  ."V 


Following  the  same  argument,  we  obtain  new  expressions  for  the  displace¬ 
ments  and  stresses.  This  is  also  presented  in  Appendix  B. 


3.2  DETERMINATION  OF  qW  AND 

The  satisfaction  of  the  stress-free  boundary  conditions  leads  to 
the  following  system  of  equations  which  can  also  be  obtained  by  applying 
the  operator  p^+k(lnp  +  B/B3)  to  (2.29)  and  setting  the  coefficients  of 
p^+k  and  p^+klnp  to  zero.  Thus  we  obtain  from  (2.29a) 

K(6)q<0)  “  0  ,  (3.8a) 

•s* 


K'(S)q«»  +  K(5)q^0)-  0  , 
and  from  (2.29b) 

K(6+k)So(k)  -  wt(a+k)St(k),  (  k  >  1  ), 

K'(6+k)g0(k)+  K(3+k)a;(k)-  Wt'(5+k)gt(k)+  Wt(«+k)at(k)  >  • 


(3.8b) 


(3.8c) 


(3.8d) 


where  R  and  Mr  are  defined  in  (2.30). 

Combining  (3.5)  and  (3.8)  with  (3.6),  we  notice  that  (3.5)  and 
(3.8)  are  the  same  as  (2.27)  and  (2.29)  for  k  <  rj,  respectively,  except 
<j^k^  now  assumes  the  role  of  a£k)  • 

For  k  *  t)  the  problem  reduces  to  solving  the  following  system  of 


equations 


-t?i  St^>at 1 ”  • 


K  (6+ 1?)^©  KCS+^go  77 


Equations  (3.9)  have  an  unique  solution  for  <j^  if 
clair  1979) 

di|K(«)|/dai  *  0  ,  (  i  -  n  -  ■  ), 

•*0 


where  n  and  m  are,  respectively,  the  order  and  rank  of 
It  is  rather  difficult  to  prove  or  disprove  equat 
ically  or  numerically.  Instead,  we  will  regard  (3.9) 
equations  for  q^  and 


(Dempsey  and  Sin- 

(3.10) 

R. 

ion  (3.10)  analyt- 
as  a  system  of  8 


,  and  solve  the  system  numerically. 


mnm 
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Chapter  IV 

APPLICATIONS  TO  COMPOSITES 

4. 1  SINGULAR  POINT  IN  COMPOSITE  MATERIALS 

We  now  consider  the  axisyaaetric  coaposite  whose  cross  section  is 
shown  in  Fig. 2.  The  two  aaterials  with  axisyaaetric  interface  SQ,  RP  are 
assuaed  to  be  transversely  isotropic  with  the  z-axis  being  the  axis  of 
syaaetry.  The  interface  aakes  an  angle  with  the  z  “  0  plane.  The 
region  SMRN  is  void  and  y  and  are  the  angles  the  two  free  surfaces 
RM  and  RN  Bake  with  the  z  “  0  plane. 

Since  equations  (2.9)  is  applicable  to  each  aaterial,  we  will  use 
the  subscript  1  or  2  separated  by  a  eonaa  to  identify  the  quantity  which 
is  associated  with  aaterial  1  or  2.  Froa  equation  (2.27)  we  notice  that 
the  undeterained  constants  for  aaterial  1  are  A^)  ...  Dq^i  w^ile  that 
for  aaterial  2  are  A^^  ...  *  The  eigenvalue  &  is  the  saae  for 

both  aaterials. 

4.2  DETERMINATION  OF  6  AND  A^  ... 

Using  (2.20b),  the  traction-free  boundary  conditions  at  angles  t  m 

and  are 

N(t.)c  0  ,  (s-1,2)  ,  (4.1) 

«.  •  «.» *  «. 

where 

T  sin^  -co tif'  0  ]  (4.2) 

JJW  0  sin^  -cos^  J  • 


The  interface  continuity  condition!  at  angle  ^  are 


NWOo  -  N(*.)o  -  0 

•W  J  M 


u  -  u  -  0  , 

M  9  1  <V  »  «  A* 


(A.  4) 


(A.  5) 


where 


u  2  , 

k-0  * 


(■-1,2)  . 


(A.  6) 


Subatitution  of  equations  (2.17)  and  (2.18)  into  equations  (4.1),  (A. A) 
and  (A. 5)  yields  the  following  system  of  recurrent  equations 


K(*)  qo(0)-  0  , 


K(5+k)3o(k)-  -t2lwt(a+k)3t(k),  (k>l)  , 


(A. 7a) 


(A. 7b) 


where 


(A. 8) 


(t  -  0,1,2  ...)  , 


in  which  the  elements  of  3^k^  are  A^k^  ...  D^k^  (s  -  1,2),  and 


K(S)  - 


0 

-av,)s0, ,u,.» 

Ol 

(A.  9) 


wt(«)- 


»«■,>  t£t , , «■, , , «■,  »*Et-2. , <*,  »J 

o 


(4.10) 


-j»,>  t£t , , »,  • « *it-i , , », .  ,  (*,. . ») 

2»,>  ts»  *£t-2. ,  <*>  •  «> 

In  (4.9)  and  (4.10),  Stjg,  Ttjg,  Utjg  are  defined  in  equations  (2.25) 
and  (>£  and  Vt  are  given  by 


G tGM)  -  (5-t+2) 


1  1  1  1  ' 

—  -  _  QtGM+D  , 

>,p,  «2P2  «,p,  «2p2 


’  1  1  1  l 

vtGM)  -  (t+1) 

0  0  0  0 


Qt(^,5+1)  , 


(4.11a) 


(4.11b) 


where  Q^^.i)  is  the  sane  as  the  one  in  (2.26). 

It  should  be  pointed  out  that  the  recurrent  relation  for  q £*)  given 
in  (2.27)  applies  to  in  (4.8).  Therefore  the  problem  reduces  to 
the  determination  of  5  and  ...  (s  -  1,2).  For  nontrivial 


solution  of  we  must  have 

I  K  (5)  |  -  0  . 


(4.12) 


This  provides  the  eigenvalue  6.  For  each  eigenvalue  6,  (4.7a)  gives  q£l 

while  (2.27)  and  (4.7b)  furnish  q^  (k  *  1,2  ...).  The  discussion  c 

^  w 

unigueness  of  solutions  for  presented  in  section  2.3  applies  here. 


4.3  DEGENERATE  CASE 

When  p,  •  p}  holds  for  one  of  the  two  aaterials  or  for  both  materi- 
als,  we  use  the  expressions  for  displaceaents  and  stresses  in  (2.38)  and 
(2.39).  Once  again,  the  stress-free  boundary  conditions  and  the  inter¬ 
face  continuity  conditions  yield  the  systea  of  recurrent  equations 
(4.7).  Equations  (4.9)-(4.10)  hold  in  which  St^g,  Tt>g,  Ut  g  are  defined 
in  equations  (2.40)  if  aaterial  s  (s  ■  1,2)  is  degenerate,  while  (4.11) 
is  replaced  by 


'  i  (a-t+Ur1*!^  *  *  ]  (4.14a) 

GtGM)  -  (a-t+2)  ,  Qt  GM+1), 

L  Pi  0i  (6-t+l)  f~*sin^'-3+4i'  *  *  J  ~ 


vt0M)  -  (t+l)  1 

~  [  0 


1  (a-t+l)f1*in^  *  * 


*  * 

.  .  Qt0A,a+l). 

0  0 


(4.14b) 


0[(^.5)  is  obtained  froa  equation  (2.26)  with  ff  •  ■  f.  As  in  equa¬ 
tion  (2.40)  the  third  and  fourth  columns  of  the  matrices  in  equations 
(4.14)  are  obtained  from  the  first  and  second  columns  by  replacing  8i 
and  f  by  -0i  and  J,  respectively. 


4.4  NUMERICAL  EXAMPLES 

We  present  two  examples  in  this  section.  In  both  exampls 
-180  ,  4,  m  90  and  ■  -90  are  taken. 
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In  the  first  exaaple  the  aatarial  1  and  2  art  both  isotropic  so 
that  3-7-  1.  We  us*  y  -  0.38,  m  -  0.3xl0*psi  for  notarial  1  and  y  - 
0.45,  m  -  0.3448xl0*psi  for  aatarial  2.  Two  negative  6'»  ara  obtainad: 


8,  -  -0.432087  ,  8}  -  -0.073520  . 


3,  is  tha  aaaa  as  p,  of  plana  strain  problaa  obtainad  by  Bogy  (1971). 
Tha  coaficiants  of  ordar  zaro  which  ara  coaplax-valuad  ara  as  follows. 
For  5,, 


a(«>  -  (0. 3446-0. 6786i)c  ,  B<o)  -  (0. 4098+0. 9122i)c  , 

0,1  10,1  1 

A<°>  -  (0.1444-0. 4991 i)c  ,  B<°>  -  (0. 4313+0. 8995i)c  , 

0t2  t  0,2  1 


and  for  iJt 


A<°>  -  (0. 1744-0. 2567i)c  ,  B<o)  -  (0. 8936-0. 0208i ) c  , 

0  9  1  2  0  ,  1  2 

A -  (0. 1576-0. 1539i)c  ,  B<o)  -  (0. 9986-0. 0533i)c  , 

0,2  20,2  2 


whara  c,  and  c2  ara  arbitrary  aultiplicativa  constants.  C's  and  D's  ara 
tha  coaplax  conjugate  of  A's  and  B's  respectively,  because  3,  and  6;  ara 


In  tha  second  axaapla  aatarial  1  is  replaced  by  a  transversely  iso¬ 
tropic  aatarial  whose  aatarial  constants  ara  (with  unit  10*psi) 

c,,  -  2.152  ,  c | j  -  0.5524  ,  c12  -  0.8115  , 

Cj,  -  34.49  ,  c44  -  0.8  . 


Tha  corresponding  p’s  and  a’s  ara 
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p,  •  0 . 1 55 1 i  ,  pj  “  1.6111  , 

а,  -  55.02  ,  ■  0.01816  . 

Again,  two  negative  5' a  art  obtained: 

б,  -  -0.484629  ,  6,  -  -0.299609  . 

Tha  coafieianta  of  ordar  zaro  ara  as  follows.  For  6,, 

A<#>  -  (0. 0486-0. 0015i)c  ,  -(-0.2669+0. 19830c  , 

0,1  10,1  ■ 

A(o)  -(-0. 0429+0. 4953i)c,  ,  B<o),  -  (0. 0351-0. 9994i)c  , 

0,2  I  1 

and  for  SJt 

A(o)  -  (0.0405+0. 03990  c  ,  B<o),  -  (0. 0330-0. 55110c  , 

A(#)  -  (0.2865-0. 2931  i)c,  ,  B<o),  -  (0. 9650+0. 2621 i) c  . 

Sinca  i+k  whara  k  is  an  arbitrary  intagar  is  not  a  root  of  (4.12)  for 
both  casas,  tha  solutions  ara  unique  up  to  tha  arbitrary  constants  c, 
and  c2. 

Tha  strass  distribution  obtainad  fro«  the  first  term  of  the  eigen¬ 
functions  associated  with  and  S2  are  plotted  in  Fig. 3  -  Fig. 6.  We 
noraalize  the  stress  by  deviding  by  the  singular  factor  pS  and  a  aulti- 
plicative  constant  c  to  sake  the  maxiaun  stress  equal  to  1. 
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Chapter  V 


EIGENFUNCTIONS  AT  AN  INTERFACE  CRACK  WITH  A  CONTACT  ZONE 

5.1  UNREALISTIC  PHENOMENON 

Ha  hava  discussed  in  Chaptar  IV  tha  strata  singularities  at  a  sin¬ 
gular  point  of  an  axisymmetric  composite  in  which  tha  fraa  surfaces  and 
tha  interface  surface  intersect.  Whan  tha  two  fraa  surfaces  make  tha 
saaa  angle  with  tha  plane  z  ■  0,  tha  free  surfaces  fora  an  interface 
crack.  Fig. 7.  Whan  the  singularity  6  is  a  complex  number,  an  oscillatory 
phenomenon  in  displacement  near  the  crack  tip  occurs  and  the  two  free 
surfaces  inter-penetrate  each  other.  To  avoid  the  unrealistic  phenom¬ 
enon,  we  assume  that  a  contact  zone  is  presented  near  the  crack  tip.  In 
a  real  composite,  the  crack  surfaces  near  the  interface  crack  tip  may, 
under  an  external  load,  open  or  close  with  or  without  friction.  The 
associated  problem  for  isotropic  composites  was  studied  for  frictionless 
contact  and  for  contact  with  friction  by  Comninou  (1977a  and  1977B) . 
Hang  (1983)  studied  the  partially  closed  interface  crack  for  anisotropic 
materials  but  the  contact  region  is  assumed  to  be  frictionless.  He  will 
use  the  asymptotic  solutions  (2. 17) , (2. 18) ,  (2.38)  and  (2.39)  to  study 
the  stress  singularities  at  the  both  ends  of  the  contact  zone. 

5.2  SINGULARITIES  AT  ENDS  OF  CCYTACT  ZONE  IN  INTERFACE  CRACK 

In  Fig. 7,  AB  is  the  contact  zone,  AC  and  AD  are  free  surfaces,  and 
BE  is  the  interface.  He  will  call  the  singularity  analyses  around  point 


■  H  I»*  <.»  ».«  M  «»  >«»  >»«  l».  i'«  *it.M  *.l  t.i 


A  and  point  B,  respectively,  Case  A  and  Case  B.  It  should  be  noted  that 
there  is  only  one  independent  angle  in  case  A  whereas  there  are  two  in 
case  B.  Let  6  be  the  angle  of  orientation  of  the  crack.  In  case  A,  \p3 
m  8,  ■  6  +  it  and  ^2  ■  6  -  it.  In  case  B,  ^3  is  arbitrary. 

Using  the  saae  notations  we  have  used  before,  we  have  the  following 
boundary  conditions  for  Case  A: 


N(*s)a  sm  0  ,  (s-1,2). 


NtyO  a  ~  Ntyr  )ff  -  0  , 

J(f,)u  ,  -  J(^3)u  +  H(^  )(7  -  0 

a*  3  ***  »  1  %  3  M>f  «  3  •  •  % 


where 


r  ■** 
m  L  0 


~ain\p  cos^r 


h(*)T  - 


s in^r tcos^+r  (sgn  t#)sin^] 
-cos2^-t (sgn  ts)sin2^r 
-coasts in^-r (sgn  tJ)cos^] 


(5.1a) 


(5.1b) 


(5.1c) 


(5.2a) 


(5.2b) 


In  (5.2b)  r  is  the  coefficient  of  friction  and  sgn  t#  stands  for  the 
sign  of  shear  traction  t#. 

Substitution  of  (2.17)  and  (2.18)  or  (2.38)  and  (2.39)  into  (5.1) 
yields  a  system  of  recurrent  equations  similar  to  (A. 7)  in  which  K  and 


W(-  have  the  expressions: 


K (6)  - 


N(f)S0  (\£  ,5) 


(5.3a) 


-N(^3)S0,2(^3,5) 


i^3)So.,^3*5)+«^3>S0.1(^,,<5)  -£(^3)G0.2(^3,5)  ’ 


0 


Wt(6)« 


N(^,)Et  (*,«) 

N«,)Et,,W„i) 

JW,)rt>1W,.«*HW,)EttI»,,« 

0 


»W,)Et  ,«  ,J) 
■*,  *  |  «  « 


In  (5.3b), 


’  Gtti(iM)+vt.lti(iM), 


(5.4.) 


(5.4b) 


in  which  Stff,  Ttfl,  Ut>#,  Gtf#  and  Vt>1  are  defined  by  (2.25)  and 
(4.11)  if  the  material  a  if  tranaveraely  isotropic  and  by  (2.40)  and 
(4.14)  if  material  s  is  a  degenerated  material. 

For  Case  B  we  have  the  following  boundary  conditions: 


N(*  )o  -  N(^  )cr  -  0  , 

mm  1  mm  9  1  Mr  *  mm  %  «  Mr 

J(^.)u  -  J(^.)u  .  +  Hty.)^  .  “  0  , 

w  1  1  mm  2  mm *  *  ■  Me »  »  M 

N(*,)<t  -  N(*,)a  -  0  , 

u  .  -  u  -  0  , 

Me  t  1  Met*  Me 


(5.5a) 

(5.5b) 

(5.5c) 

(5.5d) 


in  which  J  and  H  are  defined  in  (5.2).  Again,  (5.5)  yield  (4.7)  in 

Mr  Me 

which 


K($)  - 


N(*)S0  ,(*.«) 

N(*a)S0>1ty,,6) 

So,,^3.5> 

|>1(^,5)+H(^|)S0>1(^1,5) 

(5.6a) 
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SINGULARITIES  AT  AN  INTERFACE  CRACK  WITH  A  CONTACT  ZONE 

6.1  ST ROB 'S  FORMALISM 

An  alternative  foraulation  for  the  order  of  singularities  6  at  an 
interface  crack  with  a  contact  zone  will  be  derived  in  this  chapter.  The 
derivation  is  based  on  the  Stroh  formalism  (Stroh  1958  and  1962).  The 
Stroh  formalism,  which  has  its  origin  by  Eshelby  (1953),  provides  an 
elegant  and  powerful  method  of  treating  a  certain  class  of  two- 
dimensional  anisotropic  elasticity  problems.  Unlike  the  two-dimensional 
anisotropic  solutions  developed  by  Green  and  Zerna  (1954)  which  are 
restricted  to  plane  strain  deformations ,  the  Stroh  formalism  applies  to 
a  wide  variety  of  two-dimensional  problems  in  which  all  three  displace¬ 
ment  components  are  non-zero.  Also,  unlike  the  widely  used  Lekhnitskii ' s 
approach  (Lekhnitskii  1981)  which  breaks  down  for  orthotropic  materials 
(Ting  and  Chou  1981A)  and  requires  a  special  treatment  (Ting  and  Chou 
1981B) ,  the  Stroh  formalism  has  no  limitations  except  possibly  for  the 
degenerate  materials  in  which  the  eigenvalues  of  the  elasticity  con¬ 
stants  have  a  repeated  root  such  as  in  isotropic  materials.  The  problem 
with  degenerate  materials,  for  which  other  formalism  also  have,  can  be 
treated  separately  (Ting  1982).  However,  the  Stroh  formalism  has  since 
been  perfected  by  Barnett  and  Lothe  (1973  and  1975). 

It  can  be  seen  from  Equations  of  (2.8)  and  (2.9)  that  the  first 
order  solution  of  axisymmetric  deformation  is  the  same  as  that  of  plane- 


strain  problem.  Therefore,  the  order  of  stress  singularities  6  of  axi- 
symmetric  deformation  at  an  interface  crack  with  contact  can  be  obtained 


by  Stroh  formalism. 


6.2  BASIC  EQUATIONS 

In  a  fixed  rectangular  coordinate  system  (x1tx2,x3),  let  the 
stress-strain  law  of  an  anisotropic  elastic  material  be  given  by 

°ij  “  cijknuk,n  (6.1) 

where  repeated  indices  imply  summation,  aiy  uk,  c  ^  j  kn  are,  respective¬ 
ly,  the  stress,  displacement  and  elastic  constants  and  a  comma  stands 
for  partial  differentiation.  The  equation  of  equilibrium  are 


*ij,j  “  0  .  (6.2) 

For  the  purpose  of  the  present  analysis,  we  assume  that 

uk  -  ak26+1/(6+l),  (6.3a) 

Z  -  x,  +  px2  ,  (6.3b) 

in  which  p,  6  and  ak  are  constants  to  be  determined.  Substituting  (6.3) 
into  (6.1)  and  (6.2)  yields 

CTij  *  (cijkl  +  Pci jk2  )*kz5> 

{Q  +  p(R  +  RT)  +  p *T}  a  -  0  ,  (6-5> 


where  the  matrices  Q,  R  and  T  are  given  by 

Qik“  ci lkl  •  Rik  ’  cilk2  •  Tik  “ci2k2  • 


(6.6) 
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The  superscript  T  stands  for  the  transpose.  Equation  (6.5)  provides 
three  pairs  of  complex  conjugates  for  the  eigenvalue  p  and  the  associat¬ 
ed  eigenvector  a.  If  pa>  a^  (ffl  “  1,2,...  ,6)  are  the  eigenvalues  and 
eigenvectors,  we  will  let 

Pjb+3  "  Pm  *  Jm+3  ”  *m  »  (6.7) 


where  an  overbar  denote  the  complex  conjugate.  The  general  solution  for 
u  as  given  by  (6.3)  can  be  written  as 


D“  X 


(6.8a) 


in  which  qm  and  hj,  are  arbitrary  constants  and 

2^  -  x,  +  pmx2  -  r(cos0  +  pBsin0)  .  (6.8b) 

In  (6.8b),  r  and  8  are  the  polar  coordinates. 

Let  t £  be  the  surface  traction  on  a  radial  plane  which  makes  an 
angle  8  with  the  x,  axis.  We  then  have 

ti  “  _<7£isin0  +  <7£2COsfl  ,  (6.9) 

or,  using  (6. A),  (6.5),  (6.6)  and  (6.8b), 


(6.10a) 


where 


The  general  solution  for  the  surface  traction  can  be  written  as  (Ting 
1986) 


i  ■ ;  *  wi->. 

i 

To  derive  the  order  of  stress  singularities, 
sions  of  u  and  t  at  6  »  4>  and  (6  ±  *0  where  d  is  a 

•w  Am 

1  ing  that  Z,  of  (6.8b)  for  8-6  and  (6  t  *■)  are 

Chou  1981a) , 

I  2^(6  ±  *)  - 

and  writing  2^(6)  as 

^  Zjfd)  ■  r(ai6)  ,  ta(6)  *  cos*  ♦  pasin*  , 

(6.11)  for  6  -  6  and  (6  t  »)  bocoae 

tW  -J  r6{  ♦  h.b.fi+1(^)l  , 

3 

t(*  ±  »)  -2  ral«±i(a+1)irq_b_?£+1(*)  ♦  ,+  i(^D' 

a-1  - 

Siailar  equations  can  be  written  for  u  (*)  and  u(* 
ducing  the  new  coefficients 

qa  "  qBf5+1(d>)  ,  (a  not  sunaed)  , 

ha  *  h^f  ^  (©)  ,  (a  not  suaaed)  , 

and  noticing  that  -  -gtii*,  we  have 


(6.11) 

we  need  the  expres- 
fixed  angle.  Notic- 
related  by  (Ting  and 

(6.12) 

(6.13) 

(6.14a) 

-  -  ,  (6.14b) 

rh.bB^  <*>}  . 

t  it)  of  (6.8).  Intro- 

(6.15a) 

(6.15b) 


I 

V 
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We  a«y  consider  following  boundary  conditions: 
t(3v/2)  -  0  and  t'(-v/2)  -  0  , 


(6.19) 


t (v/2)  -  t'(v/2)  , 


(6.20a) 


u,(v/2)  -  u>/2)  , 

t,-(sgn  and  t3-(sgn  t,)kJt,  ,  at  B  -  ir/2  , 


(6.20b) 

(6.20c) 


where  k,  and  k,  are  the  coefficients  of  friction  at  and  x#  direction 
and  (sgn  t })  and  (sgn  t3)  stand  for  the  sign  of  t}  and  tJt  respectively. 
By  introducing  the  Mtrices 


J  “ 


■  1 

0 

0 

0 

0 

0  ' 

0 

0 

0 

.  G  - 

*(*gn  t,)k, 

1 

0 

0 

0 

0 

“(sgn  t,)k3 

0 

1 

(6.21) 


equations  (6.20b,c)  can  be  written  in  Mtrix  notation  as 
r_1(a>l)j[u(ir/2)  -  u'(v/2)]  ♦  Gt(»/2)  -  0  . 


(6.22) 


Substituting  (6.16)  and  (6.17)  into  (6.19),  (6.20a)  and  (6.22)  and  set¬ 
ting  d  ”  v/2  lead  to 


ei^gq  ♦  «~iiwih  “  0  , 
e’i^lV  ♦  eiiirBV  -  0  . 


(6.23) 


(6.24) 


/  /  ✓  / 


Iq  Mh  *  I  q  ♦  B  h 


(6.25) 


3 


a 


3S 


(6.26) 


j[(Aq  +  Ah)  -  (A'q  +  A/h/)]  +  G[Bq  +  Bh]  ■  0  , 

*****  ^  ^  **  **  **  *M^a  ** 

Equations  (6.23)-(6.26)  consist  of  four  homogeneous  equations  for 
q,  h,  q'  and  h'.  For  a  non-trivial  solution  the  determinant  of  the 

aw 

coefficient  matrix  must  vanish.  This  provides  the  roots  for  5.  Instead 

of  finding  the  determinant,  we  eliminate  Bh  and  B'q'  from  (6.23),  (6.24) 

*****  ***  «€ 

and  (6.25)  to  obtain 

(e1  Sir  -  e"i5,r)(Bq  -  iV)  -  0  .  (6-27) 

*****  **  **  ** 

Hence  either 

(ei6ir  _  e~iiir^  ■  2i  sin6x  »  0  ,  (6.28a) 


which  leads  to  integer  &  or 

Bq  -  B'h'  . 

"**  ** 


(6.28t) 


For  the  latter  we  substitute  h,  h'  and  q'  obtained  from  (6.23),  (6.24) 
and  (6.28b)  into  (6.26).  We  then  have 

{j  [e~^ir(AB~1-  A"i/_1)  -  ei5ir(AB_1-  aV"1)]  (6.29) 

**  «Wa  **  *****  **  ** 

-  [ei5,r  -  e"iiT]G}  Bq  -  0  . 

***  ****  ** 


It  can  be  shown  (Ting  1986)  that 

AB-1  -  a'b'"1  -  -(W  +  iD),  (6.30a) 

**  **  4W 

aI"1  -  a'b'"1  - 
**  "*» 


“(W  -  iD), 


(6.30b) 
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where  W  is  real  and  antisymmetric,  while  D  is  real,  symmetric  and  posi¬ 
tive  definite.  Equation  (6.29)  now  reduces  to 

{(JW  -  G)  -  (cot  Sir) JD}  Bq  -  0  .  (6.31) 

n>  *0*0  *0*0  *0 

For  a  non-trivial  solution  of  Bq,  we  demand  that 

0O00C 

|  (JW  -  G)  -  (cot  Sir)JD  |  -  0  .  (6.32) 

0*00*0  0*0  ***» 


If  we  expand  the  determinant,  noticing  that  W  is  antisymmetric  and  mak- 

*0 

ing  use  of  J  and  G  defined  in  (6.21),  we  obtain 


cot  Sir 


(sgn  t2)k?W,  2+(sgn  t3)k3Wia 
D, ,  +  (*8"  t  j)kJDi  2+(sgn  ts)k#Di#  * 


(6.33) 


We  see  that  when  the  friction  is  absent,  i.e.,  k2  -  k,  -  0,  S  *  -1/2  is 
the  order  of  singularity. 

If  we  apply  the  above  procedure  to  the  case  of  transversely  iso¬ 
tropic  materials  under  axisymmetric  deformation,  noticing  that  u3  ■  t,  * 
0  in  equations  (6.19)  and  (6.20),  we  obtain  the  singulari  ty 


cot  Sir 


(sgn  t2)k?W12 
D„  +  (sgn  t2)k2Di2  ' 


(6.34) 


For  isotropic  composites,  it  can  be  shown  (Ting  1986)  that 


l-2v  \~2v  1-V  1  —  V 

W12  -  -  (  - ; - )  .  D,t - -  +  -  , 

MM  MM 


(6.35a) 


w,,  “  D,,  "  D,s  "  0  , 


(6.35b) 


where  v  and  n  arc  the  Fotiion  i  ratio  and  shear  Modulus,  respectively. 
Equation  (6.33)  then  reduces  to 


cot  Sir  *  (sgn  t3)k28  , 


where  3  is  one  of  the  Dundurs  constants  (Dundurs  1970) 


m(1  -  v  )  +  m  (1  ~  *0 


(II)  Case  B 

He  way  consider  the  following  boundary  conditions; 


t(r/2)  -  t'(ir/2)  and  u(ir/2)  -  u'  (ir/2)  , 


t(3ir/2)  -  t '  (-ir/2)  , 


u,(3s-/2)  -  u'(-ir/2)  , 


t2“(sgn  tJ)kJt1  and  t3-(sgn  t3)k3t,  ,  at  6  -  3ir/2 


Equations  (6.38b,c)  can  be  written  in  Matrix  notation  as 


r'^a+DJtuO^)  -  u'  (-jt/2)  ]  +  Gt  (3*/2)  -  0  , 


(6.36a) 


(6.37) 


(6.38a) 


(6.38b) 


(6.38c) 


(6.39) 


where  k2,  kJf  sgn  t},  sgn  t3,  J  and  G  are  the  same  as  in  Case  A.  Sub¬ 
stituting  (6.16)  and  (6.17)  into  (6.37),  (6.38a)  and  (6.39)  and  setting 
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Aq  ♦  Ah  ■  A/q/+  A,h/, 


(6.41) 


eia*Bq  ♦  e“ia,rih  -  •_ia,rB/qr  +  .ia*i'h'. 


(6.42) 


j[“(e*a*Aq  ♦  e”*a*Ah)  +  (»"‘a,A,q<  ♦  e^^A^h*)], 


(6.43) 


-  G[«ia,rBq  ♦  e"ia*Bh]  -  0 


Equations  (6.40)-(6.43)  consist  of  four  hoaogeneous  aquations  for 


<g,  h,  g'  and  h'.  For  a  non-trivial  solution,  tha  determinant  of  tha 


eoafficiant  aiatrix  aiust  vanish.  This  providas  tha  roots  for  5.  Ones 


again,  by  algabraic  oparations,  wa  obtain  that  5  is  an  intagar  or  tha 


root  of  following  deterainant 


(6.44) 


JW-G  -(cot  Sir) G-JD 


-  0  , 


whara  D  and  W  ara  tha  saae  as  in  Casa  A. 


If  wa  apply  tha  shove  procadura  to  tha  casa  of  transvarsaly  iso¬ 


tropic  aatarials  undar  axisyaaetric  deformation,  noticing  that  u,  ■  t,  ■ 


0  in  aquations  (6.37)  and  (6.38),  wo  obtain  tha  singularity 


- (sgn  t ,)k  W 

cot  Sir  -  - - - 1 — *  4-  . 

D,,+(sg" 


(6.45) 


Whan  tha  friction  is  absant,  i.a.,  k?  ■  0,  5  ■  -1/2  is  the  order  of  sin¬ 


gularity. 


For  isotropic  coaposites  (6.44)  reduces  to 


cot  Sir  -  -  (sgn  t2)k2B  , 


(6.46) 


A-v  v.WA  /.vj-aV  v 


N  *.  ■.  ••  •- 
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where  3  is  one  of  the  Dundurs  constants  given  in  (6.36b).  Equation 
(6.46)  agrees  with  the  result  obtained  by  Comninou  (1978B) . 

We  have  verified  that  6  obtained  from  (6.34)  and  (6.45)  agreed  with 
that  obtained  from  (4.12)  for  Case  A  and  B,  respectively.  Since  2  is 
real,  the  unrealistic  inter-penetration  of  the  crack  surfaces  does  not 
exist. 


Chapter  VII 


CONCLUDING  REMARKS 

The  problem  of  stress  singularity  in  a  three-dimensional  elastic 
solid  that  contains  axisymmetric  notches  or  cracks  and  subjected  to  an 
axisymmetric  deformation  has  been  reduced  to  a  mathematically  two- 
dimensional  problem.  In  this  case,  it  has  been  shown  that  the  eigen¬ 
functions  for  the  singularity  associated  with  an  eigenvalue  £  contain 
not  only  the  term  p^f(^,£),  but  also  the  terms  p^+*f , (^, 5) ,  p^+^fj(^,£) 
...  where  (p,^)  is  the  polar  coordinate  with  origin  at  the  apex  of 
notches  or  cracks.  In  the  case  of  interface  crack  with  a  contact  zone, 
it  can  be  seen  from  (6.34)  and  (6.45)  that  if  £  is  an  eigenvalue,  so  is 
£+k  where  k  is  an  integer.  For  the  high  order  terms  of  the  eigenfunc¬ 
tion  solution,  equations  (4.7b)  must  be  solved.  Numerical  calculation 
showes  that  (4.7b)  has  no  solution  for  k  ■  1.  To  obtain  the  high  order 
terms  for  k  £  1  the  modified  solution  of  (3.8),  which  is  obtained  by 
dif f eretiating  (2.27)  with  respect  to  £,  has  been  used.  A  solution  for 
term  k  ■  1  is  thus  obtained  but  (3.8c)  and  (3.8d)  have  no  numerical 
solution  for  k  ■  2.  To  obtain  the  terms  associate  with  k  >  2  one  has  to 
find  the  new  solutions  by  taking  second  or  higher  derivatives  with 
respect  to  £. 

From  the  numerical  computations  we  present  the  following  conclu- 


Appendix  A 

DERIVATION  OF  EQUATIONS  (2.13) 


Substitution  of  (2.9a)  into  (2.8)  and  noticing  that  is  of 

order  p^+k+2f  we  see  that  (2.8)  is  satisfied  if 


“  0  . 


aV°)  i  &*$(<>) 

3x*  ~  P  Bz* 

a*4»(k)  l  &*$00  k-1 

-  2  x*  — - 

■■o 


,  k>l  . 


(Al) 


(A2) 


dx*  p*  J>2*  «**o  Bx 

Thus  each  term  in  (A2)  is  of  order  p^+k.  Using  (2.9b)  in  (A2)  we  obtain 


k-1 


t50At+)lt(t+l)xt*1za+k't+1+  triAt(k)2txt‘1Za+k~t+1(5+k-t+2) 

-  J,  +  J;:  , 


(A3) 


where 


J .  -  k2.xk_1"BVAf(?)1(s+l)x»Za+-‘*+11 


xrl 


k-1 


J.  -  2  xk~1-B  2  A«m) X*za+1B~,+1  (6+m-s+2) 

*  m-1 


(A4) 


(A5) 


s 

S“0  5 


By  letting  s  *  t+m-k  and  interchanging  the  summations,  we  have 


k-1  k-1 


.t-M+k-t+i 


k-1  k-1 


J,  -  t51{B-f-tAt-k-hB}xt‘1Za+k't+1(«+k-t+2). 


(A6) 


(A7) 


tea 


By  setting  the  coefficients  of  the  sane  terns  in  (A3)  to  zero,  we  obtain 


2kAkW  .  *<•> ,  <A8) 

t(t^l)A$J+  2t(a+k-t+2)A^k^  (a9) 

Vk-ti(t"k+B+1)At-k+»+1+(5+k't+2)At"k+n1- 
Equation  (A9)  can  be  rewritten  as,  by  letting  n  ■  s+k-t-1, 

t(t+l)A^+  2t(5+k-t+2)A^  (A10) 

-  X,  {sA?+k~t_1+ (5+k-t+2) A?t ,  for  t  S  k-1  . 

To  express  (A8)  in  the  forn  of  (2.13a),  we  replace  k  by  k-1  in  (A8)  to 

obtair. 

20,-l)*k<V>-  (A11) 

and  subtract  (All)  fron  (A8) .  He  have 

2kA^k)~  2(k-l)A^ki1)-  A^j15,  (A12) 

which  is  idntical  to  (2.13a).  Similarly,  we  replace  k  and  t  by  k-1  and 
t-1,  respectively,  in  (AlO)  which  reduces  to 

(t~l)tAfk_1)+  2 (t-1) (5+k-t+2)A^ki1}  (A13) 

-  X  {»A!+k_t“1+(6+k-t+2)A!!i“t_1)  . 

s-1  *  *  1 


Equation  (2.13b)  is  obtained  when  (A13)  is  subtracted  from  (AlO). 
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Jo1  At  (c,,+c1,»,pJ)(«+k-t+2)(«+k-t+l) 

-At%(|[)2(c11-c|J)  (t+1)  (6+k-t+l) 

-At\ty  (t+2)  (t+1) 

+At(k)  (c12+cia«,pp  (25+2k-2t+3) 

~A^+^2(c, ,-c, 2) (t+1)]  x*Z***-t 
+  t  A^  (c,  j+c,  3«,P?)  (6+k-t+2)  (6+k-t+l) 

-At(^2(c,  ,-c,  2)  (t+1)  (5+k-t+l) 

-At^(t+2)(t+l)]  xtzf^-tlnZ,}  +  ... 

<r<k)-  J0{[-At(k)c44(l+ni)  (i+k-t+2)  (6+k-t+l) 

-At(k^c44(l+TB,)  (26+2k-2t+3)]  *t2i+k-t 
[~A^k^  c4  4  (1+bi,  )  (6+k~t+2)  (6+k-t+l))  xtz5+k-tlnZ1}  +... 

cr^k)-t2o{[At(k^5+k-t+2)  (5+k-t+l)  +A^+k^  (t+1)  (6+k-t+l) 

+At(k)[26+2k-2t+3)+At(^(t  +  l)]  c44(l+«t)Pl  xtzf+k~t 
[At(k^5+k-t+2)+At(^(t+l)]c44(l+»,)p1  (6+k-t  +  l)xtz{+k"tlnZ,}  +  . . . 


For  the  degenerate  nateriali,  we  obtain  from  equations  (2.37),  (2.38) 

and  (2.39) 


*00-  !_  { 

+  t2o[Bt(k)zxtZ5  +  k_t+1+Dt(k)zxtZ5+k"t  +  1]  (5+k-t  +  2)}  , 


(B8) 
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u<k)«  t2o([At(kb+k-t+2)+Aji5)  (t+l)+At(k)]xtZ5“'k"t  +  1 

+  [Bj  (k}a+k-t+2)  (a+k-t+D+B^Jh+l)  (6+k-t+l) 

+B£ktea+2k-2t+3)+B^}(t+l)3zxtZ6+k't 

+  [A^  a+k-t+2) +A$}  (t+l)3  xt2a+k-t+lln 2 

+  [B^a+k-t+2)+B^(t+l)]  (a+k-t+l)zxtZa+k_tlnZ}  +  ... 

u(k)«  2  {[Af(k)3i(a+k-t+2)-Bt(k)  (3-40  (a+k-t+2) 
z  t-o  c  L 

+A^k)0i-B^k)  (3-<*v)]  xtZa+k-t+l 
+  [Bt(k)0i(«+k-t+2)  (a+k-t+1) 

+B{k)0i(2a+2k-2t+3)]zxtZa+k“t 
+  [A^3i-Bt<k)  (3-4 v)]  (5+k-t+2)xtZ5+k‘t+1lnZ 

+B^k^0i (a+k-t+2) (a+k-t+1) zxtZ^+k~^lnZl  +  ... 

-1  aik)-  2  {[(A^V-Bt^Zvei)  (a+k-t+2)  (a+k-t+1) 

2n  r  t-0  L 

+A^|[)2'y(t+1)  (a+k-t+1) +Ati^7(t+2)  (t+l)]xtZ5+k_t 
+  [Bt<k>e’(a+k-t+2)  (a+k-t+i)+B^|[)27(t+i)  (a+k-t+D 
+Bt^7(t+2)  (t+1)]  (a+k-t)zxtZa+k*t_1 
+  [(A^k)01-B^k)2v0i)  (2a+2k-2t+3)+At(^2y(t+l)]xtZ5+k't 

+  [B^k)3*((2a+2k-2t+3)  (a+k-t)  +  (a+k-t+2)  (a+k-t+D) 
+B^j[27(t+1)  (2a+2k-2t+l)+Bt(^7(t+2)  (t+1)] zxtZ5+k“t"1 
+  [(A^k)0*-B^k)2v0i)  (a+k-t+2) (a+k-t+1) 

+A^^2r(t+1)  (a+k-t+1)  +A^k^y(t +2)  (t+l)]xtZ5+k-tlnZ 

+  [B^k)31  (a+k-t+2)  (a+k-t+1)  +B^^2y(t+1)  (a+k-t  +  1) 

+B^r(t+2)  (t+1)]  (a+k-t)zxtZ5+k“t_1lnZ)  +  ... 


(B9) 


(BIO) 


(Bll) 


'V 


O 


» 


3 1 
« 


i 


j 


^k)"t?O{I(At<k)(0,'7)'Bt/(k)2,'ei)  («+k-t+2)  (6+k-t+l) 
-A^)27(t+1)  (a+k-t+l)-A^Y(t+2)  (t*l)]xtZ6+k-t 
+  [Bt(k^e*-7)  (a+k-t+2)  (5+k-t+l)-Bti^27(t+l)  (6+k-t  +  l) 
-Bt&Mt+2)  (t+i)]  (a+k-Ozxtz^-t-i 
+  tU{k>  0*-Y)-B^2vei)(2a+2k-2t+3)-A^2>(t+l)]xtZi+k-t 
+  [Bt<%*-7)  ((a+k-t+2)  (a+k-t+l)  +  (2a+2k-2t-t-3)) 

-Bt<ki2>(t+1)  (2a+2k-2t+l)-Bt^Y(t+2) 

+  [(At<k>  o*-7)-Bt(k)2v3i)  (a+k-t+2)  (a+k-t+i) 

-A^)2>(t+1)  (a+k-t+l)-A^*y(t+2)  (t+l)]xtZ,+k"tlnZ 
+  [b^b’-y)  (a+k-t+2)  (a+k-t+i)-Bt^kj27(t+i)  (a+k-t+n 

-B^Y(t+2)  (t+1)]  (a*k-t)zxtZa+k~t-1lnZj  ♦  ... 

<k) “t-o {" tAt (k)  +Bt (k) 2 (i- *-) ie"1]  (a+k-t*2)  (a>k-t*i)xtza+k_t 
-Bt(kU+k-t*2)  (a+k-t  +  1)  (a-*-k-t)zxtZa+k“t“1 
-  [A{k)+Bt(k)2(l-v)  iB_1]  (2a+2k-2t+3)  xtZa’k_t 
-Bt(k}(a+k-t*2)  (afk-t^i)^(2a-*-2k-2t^3)  (a+k-oWz4^’1-1 
-tA{k)>Bt(k)2(l-Oi8_1]  (a^k-t^-2)  (a^k-t^l)xtZa+k'tlnZ 
-B{k^a+k  -t+2)  (a+k~t  + 1)  (a+k-t)  zxtZa+k~t~  *  lnZ)  *  ... 
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-1  ffJi)-tJ0U(At/(k)8i-Bt(k)  (1-20)  (4+k-t-O) 


(B14) 


1m 


♦U^Bi-B^  (1-20)  (t*l)]  (o-k-t-Dxtz4*11'1 

♦  [B{  (k>  (4«-k-t*2)*Bji^  (t*l)3  (3+k-tO)  (a+k-t)0i«tZ6+k-t_1 
t(At(k)0i-Bt<k)  (1-20)  (24*2k-2t+3) 

♦  (At(5)Bi-  »ft)(l-20)(fl)]*tZ^k-t 
♦[Bt(k*(4*k-t*2)  (4*k-f*-l)*(24*2k-2t*3)  (4+k-t)) 

♦B^JdO)  (24*2k-2t*3)]6ir*tZ,“,k~t“1 

♦  [(At<k)Bi-Bt(k)  (l-20)(4+k-t*2) 

♦(At^Bi-Bt^(l-20)(fl)]  (4*k-fl)*tZi",k*tlnZ 

♦[B^a*k-t*2)<-B^(t*l))  (4<-k-f*-l)  (4+k-t)0iz*tZ4‘fk_t'1lnZ}  +  .  . 


5 
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Fig. 10  Cross  section  of  an  axisymmetric  composite  that  contains  a 
vertical  interface  crack  with  a  contact  2one 
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